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Abstract 

The time-scale of thermalization in holographic dual models with a chemical potential in 
diverse number of dimensions is systematically investigated using the gauge/gravity duality 
We consider a model with a thin-shell of charged dust collapsing from the boundary toward 
00 ■ the bulk interior of asymptotically anti-de Sitter (AdS) spaces. In the outer region there 
is a Reissner-Nordstrom-AdS black hole (RNAdS-BH), while in the inner region there is an 
anti-de Sitter space. We consider renormalized geodesic lengths and minimal area surfaces 
as probes of thermalization, which in the dual quantum field theory (QFT) correspond to 
two-point functions and expectation values of Wilson loops, respectively. We show how 
the behavior of these extensive probes changes for charged black holes in comparison with 
Schwarzschild-AdS black holes (AdS-BH), for different values of the black hole mass and 
charge. The full range of values of the chemical potential over temperature ratio in the 
dual QFT is investigated. In all cases, the structure of the thermalization curves shares 
similar features with those obtained from the AdS-BH. On the other hand, there is an 
important difference in comparison with the AdS-BH: the thermalization times obtained 
from the renormalized geodesic lengths and the minimal area surfaces are larger for the 
RNAdS-BH, and they increase as the black hole charge increases. 
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1 Introduction 



Quark gluon plasma (QGP) has become a subject of considerable interest after the results 
from RHIC and LHC experiments of heavy ion collisions, from where it is possible to draw 
the conclusion that the system behaves like an ideal fluid with a small shear viscosity over 
entropy density ratio (rj/s) [1, 2, 3]. The crucial outcome is that the QGP is a strongly 
coupled system. Theoretical calculations assume that the system is in thermal equilibrium, 
and it is probed at momentum scales below the equilibrium temperature T, i.e. in the 
hydrodynamical regime of the plasma. A suitable approach to the dynamical description 
of strongly coupled systems is based on the AdS/CFT correspondence [4, 5, 6]. In this 
context, many investigations have been carried out after the seminal paper by Policastro, 
Son and Starinets, where they proposed a relation between the shear viscosity of the finite- 
temperature Af = 4 supersymmetric SU (N) Yang-Mills (SYM) theory plasma, in the large 
iV limit, at strong-coupling regime, and the absorption cross section of low-energy gravitons 
by a near-extremal black three-brane [7]. The M = 4 SYM theory at finite temperature 
(just above the deconfinement transition temperature T c ) in the large N limit, qualitatively 
behaves like the strongly coupled QGP. This statement is based on a number of lattice 
QFT computations (for a recent review see [8]). The results of [7] were obtained within the 
supergravity limit. The first leading order string theory corrections to rj/s of strongly coupled 
SYM plasmas have been computed in [9]. On the other hand, electrical charge transport 
properties of strongly coupled SYM plasmas have been thoroughly studied as well. These 
include the electrical conductivity, which was firstly computed in [10] within the supergravity 
limit, while its string theory corrections were obtained in [11, 12]. The first computation of 
photoemission rates in this context was done in [13] and their corresponding leading order 
string theory corrections were obtained in [14]. 

The way in which the actual heavy ion collision process is understood involves four basic 
stages. It begins with the approach regime, where two heavy ions move toward each other 
at relativistic velocities, with kinetic energies typically about 100 GeV/nucleon. Soon after, 
once the heavy ions collide, a part of their kinetic energies transforms into intense heat and 
the plasma of quarks and gluons starts forming. This is the thermalization process. After 
that, the strongly coupled QGP is formed, and here is where the hydrodynamical description 
holds. Elliptic flow has been observed in this regime, and it has been concluded that the 
system behaves like an ideal fluid, i.e. with very low rj/s. At the end, the system expands 
and cools down, leading to the hadronization process. Since QGP in thermal equilibrium 
is strongly coupled, in principle, it is possible that the thermalization described above also 
occurs within a strongly coupled regime of QCD. Regardless of whether or not this hypothesis 
holds, it is very interesting to think about the thermalization of a strongly coupled system. 
In fact, in a series of papers the issue has been studied considering the collapse of a thin shell 
of matter, by using an AdS-Vaidya type metric [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 
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27]. Earlier studies of holographic thermalization described as a dual process of black hole 
formation have been presented in references [28, 29, 30, 31, 32, 33, 34, 35]. 

In a more realistic physical situation one should also consider the chemical potential, 
Thus, it is very interesting to study the effect of the chemical potential on the thermalization 
time scale. This can be qualitatively modeled by a collapsing thin shell of charged matter, 
leading to a thermal equilibrium configuration given by a Reissner-Nordstrom-AdS black 
hole. In this paper we carry out a detailed systematic investigation, exploring the full range of 
the chemical potential/equilibrium temperature ratio, /x/T, in diverse number of dimensions, 
for both renormalized geodesic lengths and minimal area surfaces. In a real QGP the ratio 
n/T is about 0.15 or less, thus the chemical potential effects will be limited [36, 37, 38, 39, 40]. 
For instance, at supergravity level, rj/s remains the same with finite chemical potential. 
However, it is interesting to consider the effect of the chemical potential on the thermalization 
time scale for different holographic probes, mainly to be able to investigate the behavior of 
a collapsing shell of charged matter in AdS spaces. This may have applications to other 
physical systems beyond the context of SYM plasmas, such as certain condensed matter 
systems in two and three dimensions. We would like to emphasize that although the main 
motivation for the present work has been to consider thermalization of QGPs, which concerns 
RNAdSs black holes, the systematic exploration extended to lower and higher dimensions 
presented here would perhaps motivate further studies from the AdS / CMT perspective (for 3 
and 4 dimensional backgrounds) as well as higher dimensional dual QFTs (corresponding to 
6 and 7 dimensional backgrounds). It certainly will be of interest to investigate whether these 
effective backgrounds can be lifted to ten-dimensional string theory or eleven-dimensional 
M-theory descriptions, or at least to understand the difficulties to do it. 

In section 2 we describe a model of holographic thermalization at strong coupling using 
two types of extensive geometric probes which are dual to two-point functions of QFT gauge 
invariant local operators and expectation values of rectangular Wilson loops. In section 3 
we will describe the Reissner-Nordstrom-AdS black hole in arbitrary number of dimensions. 
There is an important reason to carry out a systematic study in different dimensions, which 
is the fact that the probes of thermalization such as renormalized geodesic lengths and 
minimal area surfaces are one- and two-dimensional geometric objects, respectively. Thus, 
by increasing the number of background spacetime dimensions, it allows to explore the 
ability of these extended objects to probe additional bulk degrees of freedom. This effect 
is significant as we will show from the figures introduced in the last section of the paper. 
We will introduce a time-dependent AdS-Vaidya type metric, including charged matter in 
the collapsing shell. This is essential for the present study of thermalization, and with it 
we will study thermalization of two-point functions of gauge invariant local operators in 
the dual conformal field theory (CFT), as well as rectangular Wilson loops, in the whole 
range < fx/T < oo. This is done for QFTs in 3, 4, 5 and 6 dimensions. In the case 
of AdS 3 /CFT 2 duality there are charged extensions of the BTZ black hole [41], however 
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these are special cases that we shall not discuss here. Our results show that the effect 
of introducing chemical potential produces a delay in the thermalization time of the dual 
strongly coupled system. This effect indeed becomes more evident as the size of geometric 
probe increases. Our results are consistent with a top-down thermalization process, which 
was already observed in previous computations with no chemical potential [19, 18, 20], and is 
inherent to a strongly coupled computation. Thus, as it happens for the AdS-BH case, in the 
RNAdS-BH geometry UV modes thermalize first. This is a natural consequence of the setup 
of a collapsing thin shell, since smaller geodesies and smaller minimal area surfaces dual to 
Wilson loops thermalize before the corresponding larger probes do. It would be interesting to 
consider fat collapsing shells and see whether the top-down thermalization persists. Another 
observation from this systematic analysis is that, for sufficiently large probes, a swallow-tail 
like pattern appears in the thermalization curves, both for renormalized geodesic lengths 
and minimal area surfaces. This is a known effect in systems with two different relevant 
scales. We shall discuss about the relevance of these scales. On the other hand, this pattern 
could possibly be a consequence of a numerical artifact due to the fact that numerically 
we never have a shell with zero thickness. In addition, the swallow tails could appear from 
the breakdown to the probe approximation near the horizon, as it happens in top-down 
models when considering extremal probes instead of fundamental string probes including 
the thermalization of the string degrees of freedom [42]. After a detailed numerical analysis 
our conclusion is that the swallow tails are due to the emergence of certain solutions, from 
which only one is relevant for thermalization. 

It is interesting to note the relation between two-point functions and expectation values 
of Wilson loops with the entanglement entropy in AdS 3 and AdS 4 , respectively. Ryu and 
Takayanagi [43, 44, 45] proposed a way to calculate entanglement entropy holographically by 
considering minimal area surfaces in AdS spaces. Thus, by computing geodesies in AdS 3 and 
minimal area surfaces in AdS4, we will also be computing the evolution of the entanglement 
entropy during the thermalization of CFTs in 2 and 3 dimensions, respectively, at finite 
temperature and finite chemical potential. 

Finally, it is also interesting to note the differences in comparison with previous recent 
works. For a more complete description of investigations related to scenarios of thermaliza- 
tion of an initial field configuration far from equilibrium see [20] and references therein. In 
[19, 20] it is considered the thermalization from a collapsing shell of pressureless dust leading 
to Schwarzschild-AdS black holes in 3, 4 and 5 dimensions. The thermalization probes in 
these references are renormalized geodesic lengths, minimal area surfaces and minimal vol- 
umes. Then in [17] it is studied the evolution and scaling of the entanglement entropy for a 3 
dimensional field theory using a thermal quench, leading to an asymptotically anti-de Sitter 
Schwarzschild black hole in four dimensions. In addition, they consider an electromagnetic 
quench, which after complete thermalization, renders an extremal dyonic black hole in four 
dimensions. In this reference the entanglement entropy is studied. More recently, in ref- 
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erence [24] thermalization following a non-relativistic quench near a quantum critical point 
with non-trivial dynamical critical exponent has been studied. They consider a collapsing 
shell leading to asymptotically Lifshitz spaces. In our present work we consider all cases for 
AdSd+i spaces with d =3, 4, 5, and 6, exploring the full range of /i/T, considering different 
values of the boundary separation £. We also study the swallow tails emerging from the 
thermalization curves, for renormalized geodesic lengths and minimal area surfaces. We also 
show explicitly the shape of geodesies corresponding to zero and finite chemical potential, 
for a given boundary separation I, considering a set of different times. This allows to show 
graphically the actual shape evolution of the geodesies as the thin shell moves toward the 
bulk interior. 

2 A model of holographic thermalization at strong cou- 
pling 

In this section we describe a model of holographic thermalization of a strongly coupled CFT 
based on [19, 20]. We present a Vaidya-type metric where a thin shell of charged pressureless 
dust collapses to form a charged, asymptotically anti-de Sitter black hole. Then, we analyze 
two types of QFT operators: two-point functions of local gauge invariant operators and 
expectation values of rectangular Wilson loops, that will be used in the next sections as 
probes of holographic thermalization. 

We will generalize the Vaidya metric, which was originally introduced by P.C. Vaidya in 
1943 in order to study how the geometry of spacetime evolves in the presence of a massive 
star, which varies its mass due to the effect of radiation [46, 47]. In the present case, we 
consider d + 1-dimensional metrics that are asymptotically AdS^+i. Thus 

ds 2 = [-f(v, z)dv 2 - 2dz dv + dx 2 ] , (1) 

where f(v, z) can be an arbitrary function of v and z which approaches 1 as z goes to zero. 
The AdS radius, R, has been set to 1. The ^-coordinate is the inverse of the radial coordinate, 
z = R 2 /r, so the AdS boundary is at z = 0. In addition, x stands for coordinates at the 
boundary except the time, x = (xi, • • • , Xd-i). In particular, if we choose f(v, z) = l—m{y)z d 
we will be considering a shell (with no charge) that collapses to form a Schwarzschild-AdS 
black hole. Moreover, we can extend this geometry to include charge by defining 

f RN (z, v) = l- m(v)z d + q{vfz 2d - 2 , (2) 

where the subindex stands for Reissner-Nordstrom black hole and the expression holds for 
d > 3. 

This type of metric is a solution of the equations of motion of the action 

Srn —Vaidya — Semao-s + S ext , (3) 
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where Sema<is is the Einstein-Maxwell anti-de Sitter action 

S E MA dS = - 167r ^ (d+1) J d d+1 xV^ [1Z-F 2 + d(d - 1)] , (4) 

where the negative cosmological constant is defined as A = _^pL2 an d S ext is an external 
action, whose energy-momentum tensor and current will be specified below. 
From this action the resulting equations of motion are 

R, u - l -g, v {n-2k-F 2 )-2F, x F^ = 8nG^T^\ (5) 

-^=d p {y/=^F^) = 8nG^J7 ext) . (6) 



_ (j • '(«rt) 



The AdS-Vaidya type metric presented in Eqs.(l) and (2), together with an electromag- 
netic potential defined as 



K = --qiv^Xv, (7) 



with c = \J 2 ^-\ ' satisfy the above equations of motion provided that the external currents 



are 



8nG^T^ = z d ~ x (^L-^-m(v)-(d-l)z d - 2 



(d-l)z d 2 q(v)q(v) ) 8^5^, 



^G^J( ext) = f d l) !f 2) ^q{v)5°\ (9) 

where the dot stands for derivative with respect to coordinate v. Note that in the chargeless 
case, the energy-momentum tensor only depends on rh(v) so it is immediate to check that 
the null energy conditions are satisfied provided that m(y) is non-decreasing function of v. 
The charged case is more subtle since it depends on the specific choice of both m(v) and 
q(v). 

The functions m(y) and q(v) model the change in the mass and charge of the black hole. 
In particular, if both m(v) = M and q(y) = Q are constants, then v is the usual Eddington- 
Finkelstein coordinate, 

dz 

dv = dt ~ l-Mz« + Q^ ' (10) 

where t is the time. If we rewrite the metric in terms of t and z, then the AdS-Vaidya type 
metric corresponds to a charged black-brane type geometry, 



ds 2 = - 2 

z l 



.(1 _ Mz d + Q 2 z 2d ~ 2 )dt 2 + - ^ - H - + rfx 2 

v ^ ; 1 - M z d + Q 2 z 2d ~ 2 



(11) 
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Again, for Q = 0, we obtain a black-brane solution with an event horizon at r h = 1/zh = 



On the other hand, if m(v) = q(v) = the change of variables becomes dv = dt — dz, being 
the metric the usual AdS^+i one. Thus, if we choose m{y) = M 9(v) and q(v) = Q9(v), with 
9 being the step function, we obtain a spacetime separated into two regions by a shell: for 
v < it is an AdS^+i spacetime, otherwise it will be a Reissner-Nordstrom-AdS black hole. 
At v = we have a shell of zero thickness. Since we will carry out a numerical study it is 
convenient to consider smooth functions like 



where t>o represents a finite shell thickness. In the limit v o — > we recover the step function 
that represents a shock wave. In this limit case the null energy conditions are also satisfied, 
while in the finite thickness case the analysis has to be done with some care (see Appendix 
for a complete discussion on this). 

Once we have a metric to model the rapid injection of energy that initiates the thermal- 
ization process on the dual QFT, we have to choose a set of extended observables in the bulk 
which allow to evaluate the evolution of the system. The initial condition is that the space- 
time is AdSd + i. Then, after certain time elapses, the system reaches the thermodynamical 
equilibrium, and the bulk becomes a Reissner-Nordstrom-AdS black hole. A RNAdS-Vaidya 
type metric interpolates between these two spaces, describing the time evolution of the sys- 
tem. In order to follow this evolution we have to probe the bulk by using extended geometric 
objects, which should be computed before and after the system reaches the thermodynamical 
equilibrium. However, it is easy to imagine that not any observable can be used as a probe 
of thermalization. One may ask what happens with the expectation values of one-point 
functions of gauge invariant operators, such as the energy-momentum tensor for instance. 
The problem is that these are local operators and, therefore, cannot give information on the 
thermalization process. So, clearly one needs to consider extended probes. A possibility is 
to study correlation functions of two local gauge invariant operators on the field theory side. 

The AdS/CFT correspondence allows to understand intuitively why the expectation val- 
ues of local quantum field theory operators are insensible to the details of the thermalization. 
These observables can only account for effects due to bulk properties close to the AdS bound- 
ary, thus they do not allow to probe details from distances of the order of the thermal scale of 
the system. On the other hand, two-point functions of local gauge invariant QFT operators 
probe the bulk interior. Indeed, the AdS / CFT correspondence provides an elegant geometric 
way to compute two-point functions: under certain approximations the two-point functions 
correspond to geodesies which connect the two points where the local QFT operators are 




(12) 



(13) 



7 



inserted at the AdS boundary. These geodesies extend toward the bulk, thus allowing to 
probe a larger range of energies from the boundary field theory perspective. There are also 
other operators which are useful for the present problem: expectation values of Wilson loops 
which correspond to minimal area surfaces in the AdS bulk, as well as the entanglement 
entropy which is related to minimal volumes. 

Following [20] we first consider two-point functions of local gauge invariant QFT operators 
at constant time as probes of thermalization. In general Wightman functions are defined as 
[19, 20] 

G£(t, x; t', x') =< 0(t, x)0(t', x') > , (14) 

where O is a local gauge invariant QFT operator of conformal dimension A. In order to 
follow the evolution of these functions we are interested in the equal time (t = t') correlators. 
Then, we want to see how these functions vary at different times. When we consider strongly 
coupled QFTs these functions can only be computed analytically for dimension two. For 
higher dimensions the analytical treatment is unknown. 

Fortunately, the AdS/CFT correspondence allows to evaluate these functions when the 
operators are heavy by using geodesies in AdS spaces. As it is well-known from the AdS/CFT 
dictionary a scalar field ip(z, t, x) with mass m in d + 1 dimensions is dual to an operator O 
whose conformal dimension is A = \{d+ V d 2 + m 2 ). In general, the two-point function (14) 
in the strongly coupled regime is computed using the classical supergravity action in terms 
of ip. However, for our purpose it will be more convenient to compute it from a path integral 
as in reference [48, 20] 

< 0(t,x)0(t,x!) >= y XW AL(P) « e ~ AC > ( 15 ) 

geodesies 

where the path integral includes all possible paths connecting the points at the AdS boundary, 
i.e. (t, x) and (t, x'). In addition, L(V) is the proper length corresponding to this path. 
For space-like trajectories L(V) is imaginary. Thus, it is possible to make a saddle-point 
approximation for A ^> 1. Therefore, only trajectories with extreme geodesic lengths will 
contribute. Notice that in the last term £ indicates actual length of the geodesic between 
the points at the AdS boundary. In this way, there is a direct relation between the logarithm 
of the equal-time two-point function and the geodesic length between these two points. It is 
important to be careful while considering these approximations because the geodesic length 
diverges due to the AdS boundary contributions. Then, one can define a renormalized 
distance 5C = £ — 21n(2/^o), in terms of the cut-off z . This suppresses the divergent part 
in the pure AdS. 

The other type of non-local operators that we will be using are spatial Wilson loops, which 
are non-local gauge invariant operators in the field theory defined as the integral in a closed 
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path C of the gauge field 




(16) 



where V means ordered product, N is the number of colors of the gauge theory, and is 
the non-Abelian gauge field. Wilson loops provide information about the non-perturbative 
behavior of gauge theories, however, in general it is difficult to compute them. Using the 
AdS/CFT correspondence its computation can be done in an elegant way. The expectation 
value of a Wilson loop is related to the string theory partition function with a world sheet £ 
extended on the interior of the bulk, and ending on the closed contour C on the boundary, 



where one has to integrate over all the non-equivalent surfaces whose boundary is <9£ = C, 
at the AdS boundary. A(E) is the string action. In the strong coupling regime we can carry 
out a saddle-point approximation of the string theory partition function. In this way we can 
reduce the computation of the expectation value of a Wilson loop to determine the surface 
of minimal area of the classical world-sheet whose boundary is C. Thus, can write 



where ^4(S ) represents the area of the minimal surface. This will be a solution to the 
equations of motion of the bosonic part of the string action [52, 53]. 

These models based on the AdS/CFT correspondence allow to understand intuitively how 
the thermalization process takes place. Let us consider the sudden injection of energy in the 
dual QFT. The bulk geometry associated with this process is proposed to be described by a 
collapsing shell from the boundary toward the bulk interior. As long as the shell collapses, 
the outer region is described by a Reissner-Nordstrom-AdS black hole, while the inner region 
is still an AdS space. Now, let us use the geodesic approximation to compute the equal-time 
two-point functions. If the separation of the boundary points is small enough, then the 
geodesic cannot reach the shell at v — and, therefore, the geodesic is seen as a purely 
RNAdS-BH geodesic, i. e. for short distances in the field theory the system seems to be in 
thermal equilibrium. If we increase the separation between the insertion of the boundary 
operators, at some point, the geodesic will cross the shell, and there will be a geodesic 
refraction which will deviate the geodesic in comparison with the thermal one. Thus, we 
can understand why the thermalization proceeds from short to long distances, i.e. QFT 
ultraviolet degrees of freedom thermalize first [20]. 



Let us now consider the pure Einstein-Maxwell anti-de Sitter action (with no external fields) 
which describes a d + 1-dimensional space-time, with negative cosmological constant A = 




(17) 



< W(C) >~ e 



(18) 



3 The Reissner-Nordstrom-AdS black hole 



9 



- d ^ R P , coupled to an Abelian gauge field A M , with M — 0, • • •, d [49, 50] 

(19) 



^EMAdS = - 167r ^ d+ i / d d+1 Xy^ 9 



K-F 2 + d{d 



R 2 



This action has been introduced in Section 2, however here we recover the AdS radius R. 
The solution to the equations of motion derived from the above action leads to the Reissner- 
Nordstrom-AdS black hole metric, which can be written in static coordinates as 

dv 2 

^ 2 = -W 2 + ^+r 2 Ci, (20) 

where df^-i * s the metric on the sphere S d ~ x , and 

r 2 M O 2 
V(r) = l+jp-^ + £i- (21) 

In this metric M and Q are related to the ADM black hole mass M and the Q charge [51] 
as follows 

* = i l i^ lM - « 22 » 

Q = V2(d-l)(d-2) (^g)C, (23) 

where uia-i is the volume of the unit radius sphere S^ -1 . Also, there is a pure electrical 
gauge potential given by 

"£2+*)*' ( 24 ) 



where c = y 2 ^~^ , while $ is a constant which plays the role of the electrostatic potential in 
the region between the event horizon and the boundary of the asymptotic AdS. It is defined 
such that A t (r h ) = 0. Thus, 

< 25 > 

c ' h 

Following [49] we can consider the limit where the boundary of AdS^+i is M. d instead of 
M. xS^ 1 , the so-called infinite volume limit, which is important for the dual field theory 
discussion. This relays upon the presence of a negative cosmological constant. This limit can 
be obtained by introducing a dimensionless parameter A, and set for the radial coordinate 
r — > A 1 / d r, for time coordinate t — > A _1 / d t, for the mass and charge m — > Am and q — > 
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\(d-i)/d Qj reS pectively, and for S d - 1 we set R 2 dn 2 _ 1 ->■ A~ 2/d J^tl dx l Then > if 

we consider 

the limit A — > oo, the metric becomes 



.2 ^2 



^ = - f /(r)^ + ^- y + ^^^, (26) 

^ ' i=i 



where 



r 2 M O 2 
^) = ^-^ + ^=j- (27) 



Notice that replacing r by 1/z (setting R — 1) in Eqs.(26) and (27) we obtain Eq.(ll). 

Using the gauge/gravity duality dictionary the Hawking temperature corresponding to a 
black hole is assumed to be the equilibrium temperature of the dual QFT. For the RNAdS- 
BH metric the temperature is given by 

J-rn dr h - v 2 ; d v 3 

in terms of the black hole event horizon r\ and its charge parameter Q. This expression 
results from the usual procedure of demanding regularity at the event horizon of the Eu- 
clidean continuation of the metric. When Q = this equation leads to j3 — being the 
well-known Hawking temperature in the AdS-BH case. Then, for instance, setting d — 4 it 
reduces to the Hawking temperature of the Reissner-Nordstrom AdS 5 black hole [40], 

1 nR 2 

PRN - TP, - 2R2, • ( 29 J 

Trn r h (l - Q^- ) 

We can also consider the black hole mass parameter in terms of the radius of the event 
horizon 

M = ^ + ^2- (3°) 

We should notice that the RNAdS black hole has an extremal solution corresponding to 
T = 0. However, it nevertheless has an event horizon. This can be seen from Eq.(28), leading 
to 

r h \^={^<?#Y* • (31) 

This is a remarkable difference with respect to the Schwarzschild-AdS case, where the ex- 
tremal solution corresponds to an anti-de Sitter spacetime. Another important difference 
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is that from the AdS/CFT correspondence it is possible to understand the electromagnetic 
field in this geometry as a source of the chemical potential in the dual quantum field the- 
ory. However, the precise relation is somehow subtle since the chemical potential has energy 
units in the dual field theory ([/x] = 1/[X]) [40]. On the other hand, as defined in the 
action (19) is dimensionless, thus one has to redefine this field as = A^/R*, where i?* 
is a scale with length units. Therefore, A and /x have the same units. The only effect of 
this field redefinition is to change the Maxwell term in the action, changing the coupling as 
9ym = ^G/R*- This is useful now because in this way we can obtain the chemical potential 
in the QFT as the boundary value of the time component of the bulk gauge field. Thus, we 
obtain 



$ 1 _Q 

h 



a= lim A t = — = — — . (32) 
F ™ R* crt 2 R* 



Therefore, the black hole charge is related to the chemical potential in the gauge theory. 
depends on the particular compactification, whenever the solution could be obtained from 
string theory or M-theory. Besides the specific value of i?*, it is possible to study the whole 
range of values of /x/T given by 

t = , 4 * B2Q (33) 

Thus, if the radius of the event horizon is kept fixed, then Eqs.(31) and (32) show that it 
is possible to go from /x/T = 0, for Q — 0, to the extremal case, where /x/T — > oo, provided 
that Q satisfies Eq.(31). Recall that the mass is given by Eq.(30). 



4 Holographic thermalization with a chemical poten- 
tial 

In this section we consider how the thermalization process occurs within the model described 
here. For this purpose we will study two-point functions of local gauge invariant operators 
and expectation values of rectangular Wilson loops in the dual conformal field theory. The 
former ones correspond to geodesic lengths in the bulk theory, while the second ones cor- 
respond to minimal area surfaces in the bulk. The idea is that we start from an anti-de 
Sitter spacetime of dimension d+1. This is the holographic dual of the strongly coupled 
regime of a conformal field theory at zero temperature and zero chemical potential. Then, 
we consider a thin shell of charged dust propagating from the boundary toward the bulk in- 
terior, and collapsing, leading to a RNAdS black hole. From the dual QFT point of view this 
corresponds to a sudden injection of energy and matter, both modeled through holographic 
quenches, such that whenever the thermodynamical equilibrium is reached, the system will 
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be described within the grand canonical ensemble. As the thin shell of charged matter col- 
lapses in the bulk, it separates two regions. The outer one is a Reissner-Nordstrom black 
hole which is asymptotically AdS. The inner region is just an AdS spacetime. 

This situation is described by a dynamical metric, which is the RNAdS Vaidya type metric 
given in Eqs.(l) and (2), using the functional forms of the mass and charge of Eqs.(12) and 
(13). The conjecture is that from the holographic dual field theory this dynamical situation 
corresponds to the thermalization process of a strongly coupled system. We will see that the 
holographic thermalization is a top-down process in the sense that UV degrees of freedom 
equilibrate first. 

Now, we can actually calculate the geodesic length and the minimal area surface of the 
string attached to the boundary, as probes of thermalization for this geometry. 



4.1 Renormalized geodesic lengths 

We will be evaluating geodesic distances as function of both time and boundary separation 
length, in arbitrary number of spacetime dimensions. We will consider space-like geodesies 
between points (t,Xi) = (t , —1/2) and = (t ,£/2) in the case of AdS 3 /CFT 2 , where £ 

is the separation of the AdS boundary points. For d = 3,4, 5, 6, the orthogonal coordinates 
are fixed. For instance, for d = 4 we have (x2,x 3 ) = (x' 2 ,x' 3 ). Thus, we take as the geodesic 
parameter the first coordinate X\, which, in order to make simpler the notation, we rename 
as x. The solutions to the geodesic equations are given by the functions v(x) and z(x). 
Inserting a cut-off z close to the AdS boundary, the boundary conditions become 

z(-£/2) = z , z(£/2) = z , v(-£/2) = t , v(£/2) = t . (34) 

Also, v(x) and z(x) are symmetric under reflection x — > —x. The geodesic length is defined 

as 

c = fj=&= F ^i-M*M«) -/ w (».«M«)' , (35) 

J J-e/2 z{x) 

where /W(t>, z) was defined in Eq.(2). The prime indicates derivative with respect to x. The 
functions v(x) and z(x) minimize the geodesic length of Eq.(35), thus the problem is similar 
to the one in classical mechanics with the same Lagrangian. Since the Lagrangian does not 
depend explicitly on x, there is one conserved quantity, equivalent to the Hamiltonian of the 
system. In terms of fRN(v,z) it becomes 

n = , 1 == . (36) 

Z(x)y/1 - 2z'(x)v'(x) - f RN (v,z)v'(x) 2 

Introducing the initial conditions on the tip of the geodesic as 

z(0) = z*, v(0)=v*, v'(p) = z'(0) = 0, (37) 
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the conservation equation simplifies to 

l-2z'v'-f RN (z,v)v' 2 =^y. (38) 

The next step is to compute the equations of motion for v(x) and z(x). Though in principle 
they have a complex form, it is possible to simplify them by making use of the conservation 
equation. In particular, we can differentiate Eq.(38) with respect to x and obtain a relation 
for v"{x) and z"(x). Inserting these relations in the equations of motion, we obtain a set of 
simplified, second order differential equations for v(x) and z(x) as follows 

= 2-2v'(x) 2 f RN (v,z)-4v'(x)z'(x)-2z(x)v"(x) + z(x)v'(x) 2 d z f RN (v,z), (39) 

which is obtained from the equation of motion for z(x), while 

= z(x)v"(x)f RN (v, z) + z(x)z"(x) + z(x)z'(x)v'(x)d z f RN (v, z) 

+^z(x)v'(x) 2 d v f RN (v,z), (40) 

is obtained from the equation of motion of v(x). From this set of equations (plus the initial 
values) we can numerically obtain z(x) and v(x) for each pair of v* and z*. We can also 
obtain the equations of motion in terms of m and q, just by replacing f R ^ and its derivatives 
into Eqs.(39) and (40), leading to 

= ^dz(x) d v'(x) 2 m(v) - z(x) d v'(x) 2 m(v) - d z(x) 2d - 2 v'(x) 2 q(v) 2 

+2z(x) 2d - 2 v'(x) 2 q(v) 2 + z(x)v"(x) + 2v\x)z\x) + v\x) 2 - 1 , (41) 

= ~z(x) d v'(x) 2 m'(v) - z(x) d v"(x)m(v) - dz(x) d - 1 v'(x)z'(x)m(v) 

+z(x) 2d - 2 v\x) 2 q{v)q'{v) + z(x) 2d - 2 v"(x)q(v) 2 + 2d z(x) 2d ~ 3 v' \x)z'{x)q{v) 2 
-2z(x) 2d ~ 3 v'(x)z'(x)q(v) 2 + v"(x) + z"(x) , (42) 

where the dot stands for partial derivative with respect to v. 

The advantage of calculating the equations of motion in terms of f R N{v, z) is that they can 
be easily generalized to any other metric with that form. We can also obtain the equations of 
motion for the Schwarzschild-AdS case simply by changing f R N(v, z) for f(v, z) = l—m(v)z d , 
i.e. setting q(v) = 0. It leads to the equations for geodesic lengths presented in [20], which is 
a consistency check for our calculations and it also allows to compare the RNAdS-BH with 
the Schwarzschild-AdS one. 

At the end we are interested in evaluating the geodesic length in terms of time t and the 
boundary separation £. This information is given by 

z{£/2) = zo , v{£/2) = t . (43) 
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By using the conservation equation and reflection symmetry we can easily calculate the 
on-shell geodesic length as 



£(£ ' to) = 2 / dX 4^' (44) 

and subtract the divergent part defining 5C(£,to) = jC(£,to) — 21n(2/zo). 

Thus, we can calculate how the thermalization process occurs by considering a collapsing 
thin shell of charged dust. It is interesting to compare these results with those obtained in 
the case of a collapsing thin shell of dust with no charge, such as those described in [20]. 

In order to study the thermalization we first have to see the evolution of the geodesies as 
shown in Fig. 1. Thus, we solve the problem for a metric with an event horizon located at 
r h — l/z h — 1. The difference in comparison with the AdS-BH case is that now we have to 
set a non-vanishing charge/mass ratio. Let us consider Q/M = 1/2. This choice satisfies 
the restrictions from Eqs.(30) and (31). From Fig. 1 we can see and compare the evolution 
of geodesies with both a thin shell of charged and uncharged matter for d — 4. In both 
systems there is a similar behavior. For short times the geodesies do not change with respect 
to the pure AdS space. On the other hand, starting from t ~ 0.8, in this configuration 
with operator boundary separation £ = 2.6, differences between the charged and uncharged 
systems become apparent. Notice that within the range 0.8 < t < 1 geodesies in both 
systems change abruptly before both systems reach the thermodynamical equilibrium. Al- 
though this process is the same for the RNAdS-BH case, the actual evolution of the probes 
is different. Indeed, we notice differences which render the charged system to have a slightly 
larger thermalization time in comparison with the system at zero chemical potential. 



4.2 Minimal area surfaces 

We now focus on the thermalization of expectation values of rectangular Wilson loops. Pro- 
ceeding in a similar way as for the geodesic lengths, now we consider the computation of the 
minimal area surfaces as described before. Using the AdS-Vaidya type metric, the Nambu- 
Goto action becomes, 

a a e m R a V 1 - fRN(v,z)v' 2 -2z'v' 

A NG {to,£,R) = — / dx = . (45) 

2tt J_ e/2 z A 

We are considering boundary rectangles parametrized by the coordinates (xi,x 2 )- The 
rest of the coordinates (if there is any) at the AdS boundary are kept fixed. One assumes 
the translational invariance along x 2 , thus, we will use x\ to parametrize the functions v{x\) 
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(a) t = 0.498 - AdS-Vaidya 



(b) t = 0.498 - RNAdS Vaidya 




(c) t = 0.778 - AdS- Vaidya 



(d) to = 0.771 - RNAdS Vaidya 
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(e) t = 1-010 - AdS- Vaidya 



(f) to = 1.011 - RNAdS Vaidya 
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(g) to = 2.004 - AdS- Vaidya 



(h) to = 2.022 - RNAdS Vaidya 



Figure 1: Evolution of the shell of dust described by the AdS- Vaidya metric with = 1, compared 
with a shell of charged dust described by the RNAdS Vaidya type metric also with = 1 and 
Q/M = 1/2. In both cases the separation of the boundary field theory operator pair is i = 2.6. 
The black hole horizon is indicated by a dashed horizontal red line. The horizonal dashed blue lines 
indicate the position of the shell in each case. The time differences in each pair of figures is due to 
the fact that they have been obtained numerically, and therefore, there are minor differences. 
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and z(xi) in the AdS^+i, and we call it x. Along the x 2 direction the rectangular path on 
the boundary has length Rwl- 

As in the previous case, there is no explicit dependence on x and therefore, there is a 
conserved quantity corresponding to the Hamiltonian. The tip of the rectangular Wilson 
loop is 2*, with z'(0) = v'(0) = 0. Then, the conservation equation becomes 

l-2z'v'-f RN (v,z)v> 2 = (^)\ (46) 

The boundary conditions continue to be the same as in the geodesies case, 

z(-£/2) = z , z(£/2) = z , v(-£/2) = t , v(£/2)=t . (47) 

Next, we have to minimize the Nambu-Goto action for this geometry. Since the calcula- 
tions are similar to those for the two-point functions we will just focus on the final resulting 
expressions. The simplified equations of motion of z(x) and v(x) are respectively 

= z(x)v'(x) 2 d z f RN (v, z) - Av'{x) 2 f RN {v, z) - 2z(x)v"(x) - 8v'(x)z'{x) + 4 , (48) 

= v'(x)z'(x)d z f RN (v,z) + ^v'(x) 2 d v f RN (v,z)+v"(x)f RN (v,z) + z"(x). (49) 

Using the explicit form of f R ]\[(v, z) we obtain the following set of differential equations 

= 4z(x) d m(v)v'(x) 2 - 4z(x) 2d ~ 2 q(v) 2 v'(x) 2 - 2z(x)v"(x) 

-Sv'(x)z'(x) -Av'(x) 2 + 4, (50) 

= -z{x) d v'{x) 2 m'{v) + z(x) d m(v)v"(x) + dz{x) d ~ 1 m{v)v'{x)z'{x) 

-z(x) 2d " 2 q(v)v'(x) 2 q'(v) - z(x) 2d ~ 2 q(v) 2 v"(x) - 2dz(x) 2d ~ 3 q(v) 2 v'(x)z'(x) 
+2z(x) 2d - 3 q(v) 2 v'(x)z'(x) - v"(x) - z"(x) , (51) 

for z(x) and v(x), respectively, which we will solved numerically. These equations reduce to 
the corresponding ones presented in [20] when the charge is set to zero. 

We can again extract the information of time and boundary separation length from 

z(£/2) = z , v(£/2)=t , (52) 

and rewrite the on-shell Nambu-Goto action by making use of the conservation equation, 
obtaining 

p ri/2 2 

A(t 0: £ : R WL ) = ^ dx^. (53) 

JO Z 

Finally, we subtract the divergent part from pure AdS space by defining 

5A(t , £, R WL ) = A(t , £, Rwl) - - — . (54) 

Z 7T 
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4.3 Evolution of thermalization probes 



Having the information about the geodesic lengths and minimal area surfaces as a function 
of time, we can describe the thermalization process. In order to do it we compare 5C and 
5 A at each time with the final values 5Crn and 5Arn, obtained in the pure Reissner- 
Nordstrom-AdS black hole geometry, i.e. by setting both m(v) = M and q(v ) = Q constant. 
This procedure is repeated for different AdS-boundary separations £ and different spacetime 
dimensions. We prefer to plot the quantities 5C = 5C/£ and 5A = 5A/(Rwl£/^), such that 
the thermalization effect can be better observed, independently of £. 

Let us consider a systematic study on how the thermalization time changes within these 
geometries. If we want to have a clear comparison with thermalization in the Schwarzschild- 
AdS case without charge, it is natural to keep fixed the radius of the event horizon, say 
r h = 1. The mass of the black hole is given by Eq.(30), so we obtain 



This relation between M and Q holds from Q = up to the extremal case. From Eq.(31), 
we find 



Eq.(55) still holds for Q values larger than its extremal value. However, in that case = 1 
becomes a hidden singularity. The black hole horizon is the larger squared root of the 
function U(r) as given by Eq.(27). With this choice we can explore the whole parameter 
range of the problem. The idea is to keep fixed the event horizon radius, and consider the 
charge in the range < Q < Q ex t- Then, we will study all the range 0</i/T<ooin the 
dual field theory. 

This is shown for the two-point functions for different spacetime dimensions d = 3, 4, 5, 6 
with £ = 2, 3, 4 in Fig. 2 and for the expectation value of rectangular Wilson loops for £ = 
1, 2, 3 in Fig. 3. Firstly, we notice certain similarities in comparison with the corresponding 
Schwarzschild-AdS cases (Q = 0). In principle, the thermalization time increases as the 
separation of the CFT operators at the boundary is increased. Moreover, at a fixed value of 
the charge, the thermalization time increases as the boundary separation does. For instance, 
for £ = 2, the renormalized geodesic length practically has no significant differences for 
the thermalization time for any charge and dimension. However, interestingly, for £ = 3, 
the differences become more evident, while for £ = 4 RNAdS-BH and AdS-BH cases are 
notoriously distinct. 

One of the most interesting effects we have found is the presence of a swallow-tail pattern 
near the end of the thermalization process. This phenomenon has been reported earlier in 
[19] for thermalization of rectangular Wilson loops in AdS 4 having ji — 0, as well as for those 



M = 1 + Q 2 . 



(55) 




(56) 
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cases when the shell is quasi-stationary. The reason for the emergence of the swallow-tail 
structure is that, at certain time, there are more than just only one single geodesic minimizing 
the action. This means that the saddle-point approximation should be taken cautiously. In 
reference [19], on the other hand, it has been argued that the swallow-tail structure depends 
on the dimension of the system. However, according to our results displayed in Figs. 2 
and 3, this seems to be a universal phenomenon which depends on the penetration of the 
probes in the bulk interior. In fact, we show how for sufficiently large probes swallow tails 
are independent on the kind of probe and on the dimension of the asymptotic AdS space. 
Again, larger dimensions induce smaller variations. This effect is understood by considering 
that the fraction of the perturbed modes by the shell is smaller as the dimensionality of the 
system increases. 

5 Discussion 

In this section we discuss the numerical results obtained by using the formalism developed 
in the previous sections. We first analyze the results of renormalized geodesic lengths, then 
minimal area surfaces which are dual to expectation values of rectangular Wilson loops and 
finally we discuss about the meaning of the swallow tails in the figures. In all the figures 
corresponding to the geodesic lengths and minimal area surfaces the curves correspond to 
values of the charge from Q = to Q ext = yj with increments of 0.25. Recall that the 
mass is given by the relation M — 1 + Q 2 . The UV cut-off has been taken as z = 0.01, we 
have set r h = 1 and the shell thickness is v = 0.01. 

Figure 2 shows the results of thermalization with a chemical potential for different space- 
time dimensions indicated as AdS^+i, with d + 1 = 4, 5, 6 and 7. In all these cases we 
consider two-point functions as extensive probes. As an example we show three different 
values of the boundary separation £ =2, 3 and 4. The pictures are displayed in such a way 
that £ varies horizontally, while the spacetime dimension changes in the vertical direction. 
In every picture the vertical axis indicates the renormalized geodesic length as a function 
of t . The fluor-green curve closer to the horizontal axis shows that faster thermalization 
corresponds to the Schwarzschild-AdS black hole. The rest of curves in all pictures indicates 
increasing values of the fx/T ratio. Thus, they correspond to thermalization with increas- 
ing chemical potential described in terms of Reissner-Nordstrom AdS^+i black holes. The 
range of fx/T ratio goes from zero to infinity, and we show several curves where always the 
fluor-green ones correspond to fx/T — which coincide with the evolution of thermalization 
reported in reference [20]. As the chemical potential over temperature ratio increases the 
curves depart from the fluor-green ones, and the larger differences occur for curves corre- 
sponding to fx/T — > oo. Looking at the first column of pictures we observe that for small 
separation distances of local quantum field theory operators at the boundary, i.e. £ = 2, the 
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Figure 2: Thermalization of the renormalized geodesic lengths for a RNAdS Vaidya type metric 
with d + 1 = 4, 5, 6, 7 and the boundary separation ^ = 2,3, 4. The first curve in each picture 
(fluor green) indicates the case Q = 0, i.e. the Schwarzschild-AdS case. The last curve in each 
case corresponds to the extremal case where fj,/T —> oo. For t = 4 swallow tails appear before the 
thermalization takes place. They are shown with more detail in the corresponding insets. 
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Figure 3: Thermalizacion of the renormalized minimal area surfaces for RNAdS Vaidya type metric 
with d+ 1 = 4, 5, 6, 7, while the boundary separation is I = 1, 1.5, 2. In each picture the first curve 
(fluor green) indicates the case Q = 0, i.e. the Schwarzschild-AdS case. The rest of the curves in 

each figure correspond to Q = 0.50, 0.75, • • • , J gf^. The last case is the extremal where (jl/T — >• oo. 
Insets allow to show the swallow tails in detail. Note that in order to make easier the comparison 
among different figures the vertical and horizontal scales are kept fixed in all figures. 



21 



renormalized geodesic length becomes zero at the same t . For £ = 3 (in the second column 
of pictures) it is already possible to notice a slight enhancement of the thermalization time 
as n/T increases. This effect is much more noticeable for £ = 4, in the third column of the 
pictures. Thus, the conclusion is an enhancement of the thermalization time as the chemical 
potential over temperature ratio increases, which is more evident as the boundary separation 
increases. 

Figure 3 shows the results of thermalization with a chemical potential for different space- 
time dimensions indicated as AdS^+i, with d + 1 =4, 5, 6 and 7, for minimal area surfaces 
which turn out to be dual to expectation values of rectangular Wilson loops. Notice that we 
display values of the boundary separation, £ = 1, 1.5 and 2, which are different from the ones 
shown in Fig. 2. From the pictures we can see that i varies horizontally, while the spacetime 
dimension changes in the vertical direction as for the geodesic lengths. In the present case 
vertical axis indicates the renormalized minimal area surfaces as functions of to. The fluor- 
green curve closer to the horizontal axis shows that faster thermalization corresponds to 
the Schwarzschild-AdS black hole as before. The curves correspond to thermalization with 
increasing chemical potential described in terms of Reissner- Nordstrom AdS^+i black holes. 
The range of \xjT ratio goes from zero to infinity, and we show several curves where always 
the fluor-green ones correspond to fi/T = 0. As the chemical potential increases the curves 
depart from the fluor-green one, and the larger differences occur for curves corresponding to 
fi/T — > oo. Thermalization time is obtained when the curves reach zero. For £ = 1.5 (in 
the second column of pictures) it is already possible to notice a slight enhancement of the 
thermalization time as fx/T increases. This effect is more noticeable for £ = 2. Therefore, 
we arrive to the same conclusion as in the geodesic lengths case: there is an enhancement of 
the thermalization time as the chemical potential over temperature ratio increases, which is 
more evident as the boundary separation increases. 

Each figure 2(a-l) as well as 3(a-l) shows that the thermalization time for renormalized 
geodesic lengths and minimal area surfaces increases as fi/T increases, while the dimen- 
sion of the bulk spacetime d + 1 and the boundary separation £ of two local field theory 
operators are kept fixed. When considering the three figures for a fixed dimension in each 
horizontal sequence, there are two effects: one is that the thermalization time increases as 
the boundary separation does; the second effect is that for small £ there are no apparent 
differences in the thermalization time for any value of fi/T and, moreover, it coincides with 
the Schwarzschild-AdS black hole case. As long as one increases £ (see figures 2 and 3 (c, f, i, 
1)), the thermalization is significantly delayed as fi/T increases. This effect is understood as 
a consequence of the fact that the term proportional to the charge in fnN(v, z) has opposite 
relative sign with respect to the term proportional to the mass, so the chemical potential 
somehow produces an effect similar to the one obtained from the reduction of the black hole 
mass. In addition, it turns out that the thermalization time is slightly larger for expecta- 
tion values of rectangular Wilson loops compared with two-point functions. For the same 
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separation length in the boundary, the minimal area surface dual to the expectation values 
of Wilson loops penetrate more into the bulk, thus resulting in a more extensive probe of 
thermalization, which takes more time to measure thermal equilibrium. From the dual quan- 
tum field theory side, we expect that the system after receiving a sudden injection of energy 
and particles will be far from thermal equilibrium. We assume that its dynamics should be 
controlled by a strongly coupled Yang Mills theory. The study of the thermalization process 
cannot be done in terms of the hydrodynamic approximation, and therefore, the description 
of the process is a very hard problem. On the other hand, we can think of a rather intuitive 
argument, which although has obviously not the value of any rigorous proof, it may help to 
understand whether our results obtained using a holographic dual model of thermalization 
are qualitatively consistent with the gauge theory expectations. First notice that in order for 
the system to reach thermal equilibrium it requires both kinetic and chemical equilibration. 
The former is expected to be faster than the second one. The reason is that chemical equi- 
libration involves additional interactions. This, in principle, might delay the thermalization 
process in comparison with the situation where no chemical potential is considered. 

Also from these figures a swallow-tail pattern emerges. In the insets of the figures we 
zoom in the corresponding regions of the figures in order to see the swallow-tail structure in 
more detail. Now, let us consider the figures for a fixed £. So, increasing the AdS dimension 
we observe that the effects of chemical potential become less relevant in all cases. This is 
due to the fact that a geodesic is a one-dimensional line which is not able to account for the 
effect of other orthogonal space dimensions. Thus, if we increase d + 1, in order to observe 
the effect on the thermalization time, we have to use an extended probe having more space 
dimensions. This effect can be seen by comparison with the thermalization curves using the 
dual of expectation values of Wilson loops as extended probes. In all the figures we can 
observe that UV degrees of freedom thermalize first. 

It is interesting to analyze the cause of the swallow-tail emergence in thermalization curves 
with large boundary separation distance. In general this kind of effect is found in systems 
where there are two different scales. In the present case, there are different length scales that 
could, in principle, induce the swallow-tail pattern. For instance, the AdS radius, the inverse 
of confinement temperature T c , the UV cut-off z and the shell thickness, can all introduce 
a second scale to the problem. However, we do not find that any of these scales are relevant 
for the appearance of a swallow tail. In fact, we will show that different solutions given in 
the swallow tails are actually non-physical solutions of our problem, in the sense that they 
have not to be taken into account. This can be shown by looking at how the solutions that 
appear in the swallow tail are. As an example, we show in Fig. 4 these solutions for the 
swallow tail that appear when evaluating expectation values of rectangular Wilson loops in 
the extremal RNAdS 4 Vaidya metric. What we can appreciate is that from the three different 
solutions that generate the swallow tail at a certain time, two of them have z(x) functions 
that propagate inside the event horizon. One could nevertheless think that the appearance 
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Figure 4: Three different solutions that appear as part of the swallow tail in thermalization of 
expectation values of Wilson loops in the extremal RNAdS Vaidya geometry with d = 3. Time 
corresponds to to = 1.778 in all three cases and separation length is I = 2. The red dashed line 
corresponds to the event horizon (zh = 1), while the blue dashed one is the position of the shell. 
In the first two cases, the shell is located inside the event horizon but still have solutions inside. 
In the insets we take a closer look at the region near the horizon and consider the shell thickness 
(-0.02 < v < 0.02) in blue. 



of those solutions is due to the shell thickness that can emulate a black hole solution with a 
deeper event horizon. However, as it is shown in the insets of Figs. 4(a) and 4(b) even by 
considering the shell thickness we obtain that these solutions propagate beyond the event 
horizon. So, only the third one, that appears in Fig. 4(c), is the physical and thermalized 
solution. 

The effective model presented above is based on a double quench (for the mass and charge) 
and it has a metric which evolves in time. These features make, in principle, difficult to think 
of finding a top-down description from string theory and M-theory. However, it is certainly 
interesting to study this kind of models based on a collapsing shell of charged dust since it 
allows to study the thermalization process in a dual strongly coupled quantum field theory. 
One of the important results is that by using this geometric extended probes one finds a 
top-down thermalization regardless of the dimension of the AdS space and of whether the 
shell is charged or not. This is a consequence of the geometric setup, where for smaller 
boundary separations the dual geometric structures probe a region closer to the boundary. 
Also, as it is the case for Schwarzschild-AdS black holes there is a delay in the beginning 
of the thermalization process due to the fact that the system only feels the sudden energy 
injection at the boundary at distances of its thermal wavelength. We have carried out a 
systematic investigation by exploring the whole range of fi/T, the boundary separation and 
the dimension of the AdS spaces. Except for the case of the BTZ black hole the expressions 
require to solve them numerically. The reason to include chemical potential, using the 
grand canonical ensemble, is because realistic systems such as QGPs and condensed matter 
systems have that. In order to do so, the formalism presented in this work, studying a general 
Reissner-Nordstrom-AdS black hole Vaidya-type metric in arbitrary number of dimensions, 
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is very interesting and we hope it will motivate further studies and extensions. In all the 
cases studied the effect produced by the chemical potential is the delay of the thermalization 
process, being it more evident as the boundary separation increases. This also agrees with 
the conclusion that the thermalization of more energetic modes proceeds first in strongly 
coupled systems. 
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A On null vectors and energy conditions of the RN- 
AdS-Vaidya type metric 

In this Appendix we show how the energy conditions are satisfied for a RN-AdS-Vaidya 
type metric like the one that has been analyzed in this work as a gravitational dual to 
thermalization of field theories at finite temperature with a finite chemical potential. 

The energy-momentum tensor at each point of a manifold satisfies the so-called weak- 
energy condition [54], T MiU W fl W u > for any time-like vector W^. In particular, for any null 
vector N^, this becomes the condition T^N^N" > 0. In the case of the Vaidya metric we will 
follow [16] and consider normal null vectors to surfaces that preserve translation invariance 
in the boundary coordinates xf, that is, surfaces of constant v and z. The resulting null 
vectors are 

K = (°>|^ )> ( 57 ) 

^ = (i?V^4o) , (58) 
where we have restored units of AdS radius R. With this choice the metric and fuN becomes 

ds 2 = — [-f RN (v, z)dv 2 - 2R 2 dzdv + R 2 d* 2 ] (59) 

z 

f RN (v,z) = R 2 -m^-^^ + qiv) 2 -^^. (60) 

The normalization has been chosen such that Aq • N 2 — — 1. As the only non- vanishing 
component of the energy-momentum tensor is 

8nG d+1 T vv = {d ' R l]_ Z ^ (^R 2d - 4 m(v) ~ z d ~ 2 q{v)q{v)^ , (61) 

the null energy condition for N\ is satisfied trivially. For the energy condition depends 
on the relation between m(v) and q(v). As mentioned before, for q{y) = the condition is 
always satisfied provided that m{v) is non-negative. In the charged case, the expression can 
be analyzed in a simpler way by making use of the relation between mass, charge and event 
horizon radius given in Eq.(30). The chosen quenches satisfy 

M 

m (v) = y#W> (62) 

Q(v) = (63) 

where g(v) = 29 (v) in the zero shell thickness case and g(v) = 1 + tanh(^) in the general 
case. The null energy condition becomes 

d - 1 z d - x . , . {R Ad - & 



4 R M ~ 



10 



g(v) [t^- + Q 2 (z d h - 2 - z d - 2 g(v))j > . (64) 
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The overall function g(v) makes that the condition is satisfied if 



z < 



1 



,-d p4d-6 
h n 




d-2 



(65) 



d-2 



Q 2 



For the zero shell thickness, i.e. for the step function, the null energy condition is satisfied 
for any value of v. In the case of a smoother function, there can be a small region for large 
Q around v — 0, where the condition is not satisfied. On the other hand, we can always 
consider R and Zh to be within a parametric region where R u - Q /z 2 h d ~ 2 > Q 2 , where the null 
energy condition is satisfied for any value of v. In addition, as the limit case of zero shell 
thickness is indeed physical, we consider that it is interesting to study the whole range of 
parameters from zero charge to the extremal case. 
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